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1 Introduction 

Electroweak baryogenesis [1^ is typically described by a set of transport equations, 
which are fueled by CP- violating source terms. The source terms arise from the CP- 
violating interactions of particles in the hot plasma with the expanding bubble walls 
during a first order electroweak phase transition 0. By diffusion the sources move 
into the symmetric phase [3], where baryon number violation is fast. For walls much 
thicker than the inverse transition temperature the wall-plasma interactions can be 
treated in a WKB approximation, which corresponds to an expansion in gradients 
of the bubble profile. At first order in gradients a CP-violating shift is induced in 
the dispersion relations of particles crossing the bubble wall [3]. A (semiclassical) 
force results, different for particles and antiparticles, which creates a non-zero left- 
handed quark density in front of the bubble. The weak sphalerons partly transform 
this left-handed quark density into a baryon asymmetry. 

The WKB approach has been widely used to study electroweak baryogenesis in 
various extensions of the standard model (SM) El El El 13 GDI El ■ (An alternative 
approach was followed in ref. 521-) In the simplest manner, the WKB dispersion 
relations were computed by solving the one-particle Dirac equation to first order in 
gradients in the CP-violating bubble wall background. In a more rigorous treatment 
similar dispersion relations were also derived in the Schwinger-Keldysh formalism 

Comparing the dispersion relation of a single Dirac fermion obtained from the 
Dirac equation j7] to that of the Schwinger-Keldysh formalism ^3] , we observe that 
the CP-violating part of the latter is somewhat enhanced. In this letter we show 
that this mismatch disappears when the result gained from the Dirac equation is 
correctly boosted to a general Lorentz frame. Thus in the case of a single Dirac 
fermion, the full Schwinger-Keldysh result can be obtained in a much simpler way. 

To demonstrate the numerical significance of this effect we recompute the baryon 
asymmetry in the SM augmented by dimension-6 operators Lllij. Using the correct 
dispersion relations enhances the baryon asymmetry by a factor of up to about 2. 
We also improve on the transport equations, keeping scatterings with W bosons 
at a finite rate, which considerably reduces or enhances the baryon asymmetry, 
depending on the wall velocity. We also show that the position dependence of certain 
thermal averages in the transport equations has a substantial impact on the baryon 
asymmetry. Finally, we investigate to what extent the CP-violating source terms 
are influenced by CP-conserving perturbations in the plasma, an effect that turns 
out to be negligible. In total, depending on the model parameters, our refinements 
can increase the baryon asymmetry by a factor of up to about 5. 

2 The semiclassical force 

We consider a single Dirac fermion, such as the top quark. Its mass changes as it 
passes the bubble wall. Once the bubble has sufficiently grown, we can approximate 



the bubble wall by a planar profile. The profile is kink-shaped and characterized 
by a wall thickness L^. The problem is most simply treated in the rest frame of 
the bubble wall. In the presence of CP-violation, the fermion mass term can be 
complex, i.e. Re(A1) +i'j^Iin{A4), where 

M = m(2)e^^(^) (1) 

and z is the coordinate perpendicular to the bubble wall. 

For a particle with momentum much larger than L~^ we can solve the Dirac 
equation using a WKB ansatz 

and expand in gradients of Ai. Here Pcz is the canonical momentum along the z 
direction. To simplify the solution we have boosted to the frame where the momen- 
tum perpendicular to the wall is zero. Since the typical momentum of a particle in 
the plasma is on the order of the temperature T, this approach is valid for thick 
bubbles, i.e. TL^ ^ 1. Note that at this stage the fermion is treated as a free 
particle. Scatterings with particles in the plasma will be incorporated later on by 
means of the Boltzmann equation. 

As shown in ref. the dispersion relation is, to first order in gradients 

^ = ViPcz - acpY + m? =F — , (3) 

with 6' = dz9, acp = a' ± |-, and s = 1 (—1) for z-spin up (down). The upper 
(lower) sign corresponds to particles and antiparticles, respectively, which this way 
get different dispersion relations. The additional phase a is related to an ambiguity 
in the definition of the canonical momentum, when replacing \1/ —^ e*"^^-*\l/. It was 
the main result of refs. E] that this ambiguity disappears when all quantities are 
expressed in terms of the kinetic momentum rather than the canonical momentum. 

In ref. jl] the dispersion relation was used to compute the semiclassical force, 
which was then generalized to a Lorentz frame with finite momentum parallel to the 
wall. The point we make in this letter is that first eq. (jHI) should be boosted to the 
general frame and all further manipulations should be carried out later on. We will 
demonstrate that this way the dispersion relation of ref. ^3] is correctly reproduced. 

Since Lorentz invariance is not broken parallel to the wall, we simply have to 
replace cu^ -^ uj'^+p'^+Py. Note that parallel to the wall we do not have to distinguish 
between kinetic and canonical momentum, i.e. Pcx,y = Px,y The dispersion relation 
Q turns into 

where 



^0 = \J iPcz - acpY + pI +pI + w? 
i^oz = ViPcz - acpY + m?. (5) 



In the limit ujq = uqz we are back at the old result. In the following we show that 
when written in terms of the kinetic momentum the dependence on acp still drops. 
The physical kinetic z-momentum is given by pz = ujVgz, where Vgz, the group 
velocity of the WKB wave-packet in the z direction, is given by 

\dpczj, u;o V 2 iu^iuoz J 

The kinetic momentum then is 

Pz = (Pcz -a) (it .^) . (7) 

We can use this expression to replace the canonical momentum in the dispersion 
relation (JH). To stress the difference, we introduce a new symbol, E, to denote 
energy expressed in terms of the kinetic momentum. Defining 



-^0 = \ pI + Pl+pl + m^ 



Eoz = VpfT^, (8) 

we obtain, to first order in gradients 

E = Eo±AE = 

Notice that the ambiguity related to acp has disappeared. For the group velocity 
we now find 

From the canonical equations of motion we can compute the force acting on the 
particle 



Pz = ^Vgz 



(11) 

where we have used the fact that u is constant along the trajectory. Performing the 
partial derivatives and replacing the canonical by the kinetic momentum, we finally 
obtain 

2Eq 2EqEqz AE^Eqz 
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Thus particles and antiparticles experience a different force as they pass the bubble 
wall. This CP-violating part of the force is second order in derivatives. There is 
also a CP-conserving part, which is first order in derivatives. 

Our expressions for the dispersion relation 0, the group velocity (fTUjl . and the 
semiclassical force (I12|l agree with the results of ref. Jl], demonstrating that for a 
single Dirac fermion the full Schwinger-Keldysh result can be obtained in a much 
simpler way by means of the Dirac equation. This is the main result of this letter. 

In the special case Eq = Eq^, i-e. when the particle has no momentum parallel 
to the wall, our results agree with those of ref. [Zj. For a relativistic particle in 
the plasma Eoz contains only roughly a third of the total energy. Keeping correct 
track of the factors Eq^ enhances the CP-violating part of the dispersion relation 
and the force term by a factor of up to about 3. For non-relativistic particles 
the effect is smaller. This factor has been neglected so far in computations of the 
baryon asymmetry based on the WKB approximation of the Dirac equation. We 
will demonstrate this enhancement in a numerical example in section 4. 

In the next section we discuss the impact of the CP-violating force on the trans- 
port equations of particles in the plasma. In a chiral theory, as the SM, interactions 
are related to the chirality of a particle rather than its spin. Thus it is convenient 
to label particles in terms of helicity A, which is close to chirality for relativistic 
particles. We then have to replace the spin by s = Asign(p^) in eqs. (jHl), (|Tnjl and 

3 Transport equations 

In the derivation of the transport equations we closely follow ref. f^. A crucial 
assumption made in that work is that it is the kinetic momentum that is conserved 
in the scatterings of WKB particles. The equilibrium phase space distributions 
should therefore also be written in terms of the kinetic momentum. In the wall 
frame this reads 



where (3 = 1/T and 7^, = l/y 1 — f^, and plus (minus) refers to fermions (bosons), 
respectively. We model the perturbations from equilibrium caused by the passage 
of the bubble wall with a fluid-type ansatz 

The chemical potentials fii{z) describe a local departure from the equilibrium parti- 
cle density. The perturbations 6fi model a departure from kinetic equilibrium and 
allow the particles to move in response to the force exerted by the bubble wall. 
They do not contribute to the particle density, i.e. J d^p 6fi = 0. To second order in 
derivatives, we have to distinguish between particle and antiparticle perturbations. 



which we can expand as 

/^i = /^i,le + /^j,2o + /^«,2e; ^ fi = ^ fi,le + ^ fi,2o + ^ fi,2e- (15) 

Notice that the second order perturbations have CP-even and CP-odd parts, which 
we treat separately. 

Let us now concentrate on the Dirac fermion of the last section, so that we can 
drop the index i to simplify the notation. We expand its distribution function to 
second order in derivatives as 

/ ~ fo,v^ + fo,v,, ilw^E - file - H2o - f^2e) 

+Sfle + Sf2o + Sf2e. (16) 

Here /o,„„ denotes the equilibrium distribution (fT^ where E is replaced by Eq, and 
fovi, — {d,/dEo)fQ.u^. The dependence on the wall velocity is taken exact at this 
stage. 

The evolution of / is governed by the Boltzmann equation 

L[f]^{zd,+p,d,Jf = C[f]. (17) 

We look for a stationary solution, so that the explicit time derivative drops. Plugging 
the ansatz (J16j) into the Boltzmann equation, taking i and p^ from eqs. (jl(J|) and 
(|T^. and subtracting the results of particles and antiparticles, we obtain for the 



flow part 
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Note that the second order perturbations present differences for particles and an- 
tiparticles, i.e. fi2 = /i2o — ^20, the same as for 6/2- The CP-even parts drop. For 
the first order perturbations we take /ii = /iie + fiie, etc. 

We average the Boltzmann equation over momentum, weighting it by 1 and pz/E^. 
We also expand in the wall velocity, keeping only the linear order, i.e. fo^^^ ~ 
/o + VwPzfl)- We then obtain 

VwKifi'2 + Vu,K2{m^)' ^.2 + u'2 - (C[/]) = Sf, 






Se + Su 



(19) 



with the source terms 



Se 



1^2 , 1 



K79'm^fi[ 

-v^Ksim^ey + v^Kge'm^m'^y 

-Kiom'^e'u[. 



(20) 



The primes again denote derivatives with respect to z. The sources S'^^„ are related 
to the first order perturbations. Notice that these are first order in Vw Formally, 



S^^u are one order higher in gradients than Sq. It will turn out that they indeed con- 
tribute only a small fraction to the total source term. After momentum integration 
we normalize the resulting equations by the average of the massless Fermi-Dirac 

distribution 

/ d^p X{p) 



{X) 



(21) 



This normalization we also use for bosons to keep the interaction rates for fermions 
and bosons equal. The plasma velocity we define as 



The thermal averages read 
K, = - 
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(23) 



The averages Ki are related to averages involving 6/2- Since we do not know the 
momentum dependence of 6/2, we make the additional assumption that these aver- 
ages factorize and then use eq. (|^. e.g. {plSf2) ~ [plEofo^y^]u. We normalize these 
averages by the massive distribution of the boson or fermion under consideration, 
i.e. f d^p fo^v^. Since there is some arbitrariness in this procedure, we will test the 
impact of these averages, which turns out to be smalF. 
The collision integrals read [7] 



(C[/]) 



r 



incl 



$^/^*,2 







■c[/] 



r*°^M2, 



(24) 



^Depending on how we precisely treat the averages involving 6/2, there can also arise a source 
term of the form {rn^yO'ui. We do not discuss it in more detail since the source terms related to 
the first order perturbations are small anyway. 



where F"^'^' and r*°* are the inelastic and total interaction rates, respectively. The 
negative sign in front of r*°* is related to our sign convention for the plasma velocity 
^. 

The transport equations of the first order perturbations look very similar to 
eq. (UnD 

-Kifx[+v^K5u[ + v^ke{m^)'ui + T'°\i = 0. (25) 

The source term is now first order in derivatives and CP-even. Note that here also 
the quite large annihilation rates enter in r"'''. 

In eqs. p9|) and (j^3|) we can approximately eliminate the plasma velocity to 
obtain diffusion equations for chemical potentials. From the coefficient of the /i" 
term we can read off the diffusion constant as [7, 

D = -^. (26) 

Our source terms (j2Uj) agree with those obtained from the Schwinger-Keldysh 
formalism. However, in ref. [T3j there is one extra source term, related to the gradient 
renormalization of the Wigner function. This term seems to be missing in the Dirac 
equation approach. It is of order m^, like the i^g-part of Sg. We will demonstrate 
in the next section that these terms are subleading. 

4 Top transport: an example 

We now apply the general results ()19|) and (f^ to top transport in an effective SM 
with dimension-6 operators [121 UHl El UM ilfl^ • The model contains a single Higgs 
doublet, whose potential is stabilized by a 0^ interaction 

V(0) = -^02 ^ A 4 1 6_ 

This potential has two free parameters, the suppression scale M of the dimension-6 
operator and the quartic coupling A. The latter can be eliminated in terms of the 
physical Higgs mass mn- Since the potential is stabilized by the 0^ term, A can 
be negative. In this case a barrier in Higgs potential is present at tree-level, which 
triggers a first order electroweak phase transition. Computing the 1-loop thermal 
potential, it was shown in ref. [12] that the phase transition is strong enough to avoid 
baryon number washout, i.e. ^ = (0)t^/Tc > 1.1 |19.^, if M < 850 GeV and tuh = 115 
GeV. Taking M = 500 GeV, a strong phase transition is present for rriH < 180 
GeV. Thus the model allows for a strong phase transition in a large part of its 
parameter space. In ref. ^U] also the wall thickness has been determined, showing 
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that 3 < LyjTc < 16. The gradient expansion discussed in section 2 is therefore 
justified in almost the full parameter space. The thinnest walls correspond to a 
very strong phase transition, ^ ~ 3, where the model is close to metastability of 
the symmetric phase. In the following we will approximate the wall profile by a 
hyperbolic tangent, 0(z) = (fc/2)(l — tanh(z/Ltt,)). 

Dimension-6 operators also induce new sources of CP-violation. In addition to 
the ordinary Yukawa interaction of the top quark, yt^f^q^, we have an operator 
(a;t/M^)($^$)<l>t'^g3 [12]. We denote the relative phase between the two couplings as 
(ft = a,Yg{yx*). Then the top develops a position dependent complex phase 9t along 
the bubble wall (biz), with 



tan et{z) ^sm^t^^ 



(28) 



So all necessary ingredients are present to apply the formalism discussed in the 
previous sections. 

For the generation of the baryon asymmetry, the most important particle species 
are the left- and right-handed top quarks, and the Higgs bosons. We will show that 
the latter have only a minor impact. We ignore leptons, which are only produced 
by small Yukawa couplings. In contrast to all previous investigations we include 
the W scatterings with a finite rate Tw- This procedure allows us to study the 
perturbations of bottom and top quarks separately. The top quark source is no 
longer locked to the bottom degrees of freedom, which would lead to a larger or 
smaller baryon asymmetry, depending on the wall velocity. The other interactions 
we take into account are the top Yukawa interaction, Ty, the weak and strong 
sphalerons, T^g and Ygg, the top helicity fiips, Tm, and Higgs number violation Th- 
The latter two are only present in the broken phase. 

In a first step we compute the left-handed quark density, assuming that baryon 
number is conserved. Later on, the left-handed quark density will be converted into 
a baryon asymmetry by the weak sphalerons. The transport equations for chemical 
potentials of left-handed SU(2) doublet tops fj,t,2, left-handed SU(2) doublet bottoms 
fib,2, left-handed SU(2) singlet tops fit'=,2, Higgs bosons fih,2, and the corresponding 



plasma velocities read 

-3r,,[(l + 9iri,t)/it,2 + (1 + 9Ki^b)l^b,2 + (1 - 9iri,t)/ite,2] = 3i^7,t^>t/^;,i 

3VwKi^bf^b,2 + 3Mfe,2 
— 3ry(/i{, 2 + /^t':,2 + Ai/1.,2) — 3rvi/(/ife,2 — fJ't,2) 

-3Tss[il + 9K^,t)l^t,2 + il + 9K^,b)l^b,2 + il-9K^,t)l^t^,2] = 

3v.ujKi^tlJ''tc2 + 3vu,K2,ti'mfyfit^^2 + 3mjc 2 
-3r,,[(l + 9iri,t)^t,2 + (1 + 9Ki^b)fib,2 + (1 - 9iri,t)/itc,2] = 3X7,t^>?/^;c,i 

2VyjKi^hfJ''h,2 + 2Mft,2 

-3ry(/ii,2 + yU6,2 + 2yUfC2 + 2/i;,,2) -2r,,/i/i,2 = (29) 

- 3ir4,t/iJ 2 + 3v^k^^tu't^2 + 3v^kQ^tim'fyut,2 + 3rf*Mt,2 = 
= -3t;^ir8,t(m2eO' + 3v^K^^tO[ml{ml)' - 3kio,tmX<,t 

-37^4,6^2 + 3t^«,i^5X2 + 3rrV2 = 

-3i^4,t/ijc_2 + 3t;t„^5,tMjc_2 + 3v^KQ^t{rnfyutc^2 + 3rf*Mtc^2 = 
= -3vy,Ks,t{mf6^y + 3vy,Kg^tO't'mf{mfy - 3kio^tmf6[u[^c 

-2Ki^hli'h,2 + '^Vu,k^,hu'h^^ + 2TTuh,2 = (30) 

In eqs. (IHUj) F^y can be neglected since the plasma velocities of t and b are damped 
by the much faster gluon scatterings. We have used baryon number conservation 
to express the sphaleron interaction in terms of fit,2, l^b,2 and /i^c^ (201 • ^ possible 
source term for the bottom quark is suppressed by {rrib/mt) and therefore neglected. 
The first order perturbations of t can be computed from 

3v^Ki^tfJ''t,i + 3f^-ft'2,t("^i)Vt,i + 3Mi_i - 3rfVt,i = 3f^K3,t(m^^)' 

-3/^4X1 + 3t^«,i^5,t<i + 3^-^6,t(^?)Vi + 3rfVi = 0. (31) 

The damping of Ut^i is dominated by gluon annihilation, the rate of which we have 
approximated by r^°*. Other scatterings have been neglected. To this approximation 
the chemical potentials of t and f^ are identical. This guarantees that no direct source 
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for baryon number is induced. Such a source can be generated if fit^i 7^ /itc,i- It 
leads to spurious effects in the baryon asymmetry. Its appearance shows that an 
inconsistent approximation pattern has been used. 

We can now compute the chemical potential of left-handed quarks, fiB^ = A''<ji,2 + 
/^52,2 + (/^i,2 + /^b,2)/2. Assuming again baryon number conservation, we obtain 

1 1 

A^Bi = 2(1 + 4i^i,t)/^t,2 + 2 ^^ + 4^i.fe)/^M - 2i^i,t/iic,2. (32) 

The baryon asymmetry is then given by [7j 



Vb = — 



ub 405r 



ws 



S Ati'^VuiQ^T 



/ rf2/iBJ2)e""^ (33) 

Jo 



where is T^s the weak sphaleron rate and v = 45ru;s/(4fu,). The effective number 
of degrees of freedom in the plasma is g^, = 106.75. In eq. (jH^ the weak sphaleron 
rate has been suddenly switched off in the broken phase, 2; < 0. The exponential 
factor in the integrand accounts for the relaxation of the baryon number if the wall 
moves very slowly. Note that we have performed our computation in the wall frame. 
Therefore, strictly speaking eq. (^^ gives the baryon asymmetry in that frame. To 
first order in v^ this is identical to the baryon asymmetry in the plasma frame. 

In our numerical evaluations we use the following values for the weak sphaleron 
rate [21], the strong sphaleron rate ^^, the top Yukawa rate [20] , the top helicity 
flip rate, the Higgs number violating rate [20], the quark diffusion constant jl| and 
the Higgs diffusion constant |7j 

r^, = 1.0 X 10-^T, r,, = 4.9 X 10-^T, 



r„ = 4.2 X lo^^T, r 



y 



63T 



mw{z,T) D =- 

50T ' '^ T' 



Dh = f . (34) 

We use eq. (12 (ij) to infer the total interaction rates from the diffusion constants. 
In this procedure we evaluate the thermal averages at z = 0, i.e. in the center of 
the bubble wall. The W scatterings we approximate as T^ = T)^^. The bottom 
quark is taken as massless, and the Higgses we count as 2 massless complex degrees 
of freedom. The rates of eq. ()34j) have been computed in the plasma frame. We 
assume that, to leading order in v^, they can also be used in the wall frame. 

To demonstrate the relevance of the various contributions to the full transport 
equations, we compare the baryon asymmetry computed in different approximations 
for two typical parameter settings. We take \xt\ = 1 and maximal CP violation 
simpt = 1- Fig. H shows rjB as a function of the wall velocity Vw The other 
parameters we have chosen as ^ = 1.5, M = 6 and L„, = 8. These values correspond 
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Figure 1: The baryon asymmetry as a function of v^ for ^ = 1.5, M = 6 and Ly 
The labehng is explained in the text. 
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Figure 2: The baryon asymmetry as a function of v^ for $^ = 2.5, M = 6 and L„ = 3. 
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to a setting where the baryon asymmetry is close to the observed value tjb = (8.9 ± 
0.4) X 10-" 123 111- In Fig. 12 we use ^ = 2.5, M = 6 and L^ = 3, i.e. a very 
strong phase transition with a small wall width. In both figures the bold solid line 
(a) indicates tjb using the source terms Se and keeping the full 2;-dependence of the 
thermal averages (jSHI)- 

In (b) we drop the space-dependence of the thermal averages. We rather evaluate 
them at the center of the bubble wall, i.e. Ki^t{z) = Ki^t{z = 0). Formally, the 
space-dependence of the thermal averages is a higher order effect in gradients. But 
this approximation considerably underestimates the baryon asymmetry, especially 
for small wall velocities and thin bubble walls. The full z-dependence reduces the 
impact of the wall velocity on rjB- 

The long-dashed line (c) shows the result when we resubstitute Eq^ -^ Eq, going 
back to the dispersion relation of ref . . This would considerably reduce the baryon 
asymmetry, in particular for weaker phase transitions (Fig. ^. 

Neglecting the Higgs bosons in the transport eqs. ^I^ and (IHUj) leads to a reduc- 
tion of r]B by ~ 10% (d), almost independent of the wall velocity and the strength 
of the phase transition. 

Taking the W scatterings to equilibrium (e) has a substantial effect on the re- 
sulting baryon asymmetry, especially for strong phase transitions. In Fig. |2lit over- 
estimates rjB by a factor of almost 2 for v^ < 0.1. For large wall velocities there is 
an underestimate of t]b by a similar size. Keeping W scatterings finite results in a 
much milder f ^-dependence of the baryon asymmetry. 

The dash-dotted line (f) adds the contributions of S"^ + Su to line (a). The effect 
of these source terms is quite small, consistent with the fact that they are of higher 
order in gradients. They enhance the baryon asymmetry in the whole f ^-range only 
by a few percent. 

Line (g) shows the effect of switching off the terms proportional to K^ and Kq. If 
these terms are neglected, the final result is reduced by a contribution proportional 
to the wall velocity. It demonstrates that the precise treatment of the averages 
involving 6f has only a minor impact on the baryon asymmetry^. 

Altogether the examples demonstrate that the leading contribution to rjB comes 
from the source Sq. The baryon asymmetry gets considerably enhanced by using the 
dispersions relation with the correct factors of £^02 and keeping the space-dependence 
of the thermal averages. The finite W scattering rate has a sizable effect, the 
direction of which depends on the wall velocity. The resulting u^-dependence of 
the baryon asymmetry is rather mild. The baryon asymmetry grows slowly with 
increasing f^ and reaches a maximum at v^ — 0.2-0.3. Taking the Higgs bosons or 
the S^ + Su sources into account is less important. Their effect is not larger than 
typical uncertainties from higher order terms in the gradient expansion. 

The source Sg consists of two parts, proportional to Kg and Kg. The latter 



^Numerically there is also not much difference to the prescription used, for instance, in ref. |1U|. 
where plasma velocities were included in the fluid ansatz, rather than using a general Sf. Then, 
for example, the W2 term in eq. H19|l obtains an additional coefficient ^^ 1.1. 
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Figure 3: The baryon asymmetry in the SM with low cut-off for two different Higgs 
masses as a function of M (in units of GeV) for v^ = 0.01 (sohd) and v^ = 0.3 
(dashed). The horizontal lines indicate the error band of the observed value. 

has an additional factor m^, leading to an extra suppression, in particular for weak 
phase transitions. For instance, taking f^ = 0.1 and the parameter set of Fig.^ the 
i^g-part contributes only about 15% to the total baryon asymmetry. As indicated 
earlier, there is an extra source term in ref. f^, which is related to the gradient 
renormalization of the Wigner function. It also has an extra factor of m? and 
therefore should also be sub- leading in our case. 

Of course, the baryon asymmetry also depends on the precise values of the in- 
teraction rates (J3l|l . For instance, reducing the quark diffusion constant by 10% 
leads to an about 7% reduction in the baryon asymmetry (taking Vw = 0.1 and the 
parameter set of Figure H}. Changing Vw by 10% affects ttjb to less than 1%, even 
for Vw ~ 0.01, where the impact of the W scatterings is particularly large. 

FigureEldisplays the baryon asymmetry in the SM with a low cut-off as a function 
of the cut-off scale M. We consider two different Higgs masses ruH = 115 GeV and 
rriH = 150 GeV and two wall velocities Vi^ = 0.01 and 0.3. For each value of M the 
corresponding strength of the phase transition and bubble width are computed as 
in ref. ^U]. As expected tjb increases rapidly with decreasing cut-off scale M. The 
asymmetry has only a minor dependence on the wall velocity. In both cases it is 
possible to generate the measured baryon asymmetry for a reasonably small value 
of M. 
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5 Conclusions 

We have studied aspects of top transport in electroweak baryogenesis. We have 
computed its CP-violating source term in the WKB approximation, using the one- 
particle Dirac equation in the wall background. When the top dispersion relation is 
correctly boosted to a general Lorentz frame, the Schwinger-Keldysh result |13| IT^ 
for the semiclassical force term is obtained in eq. (fT^ . The CP-violating source 
term is enhanced with respect to ref. pj. We have only considered the case of a 
single Dirac fermion, but our results should simply generalize to mixing fermions, 
such as the charginos in the MSSM. 

In our computation we cannot obtain the extra source term of ref. ^3] , which is 
related to the gradient renormalization of the Wigner function. In the case of top 
transport this term is subleading since it is of order m'^. In our approach, of course, 
we also cannot obtain source terms related to quantum mechanical oscillations be- 
tween different fermion flavors. In the case of the top quark this effect is obviously 
not present, but it can be relevant for the charginos in the MSSM |25j . 

We have demonstrated the numerical significance of the corrected dispersion 
relations in the SM augmented by dimension-6 operators. This effect alone enhances 
the baryon asymmetry by a factor of up to about 2. We have also improved on the 
transport equations, keeping scatterings with W bosons at a finite rate. Depending 
on the wall velocity and the wall thickness, putting the W scatterings to equilibrium 
(as was done so far in the literature) can increase or decrease the baryon asymmetry 
by a factor of 2. It would be interesting to study the impact of this effect in 
supersymmetric models, where the SU(2) supergauge interactions have been put 
to equilibrium as well. 

We have shown that the position dependence of the thermal averages in the 
transport equations has a substantial impact on the baryon asymmetry, even though 
it is formally a higher order effect in the gradient expansion. Finally, the influence 
of the Higgs bosons on transport turned out to be small, as is the contribution of 
the sources S^^u (EHl)- In total, depending on the model parameters, our refinements 
can increase the baryon asymmetry by a factor of up to about 5. 

The rather large impact of the precise treatment of the W scattering rate and 
the space-dependence of the thermal averages probably indicate that there is still a 
substantial uncertainty related to transport. 

In a forthcoming publication we will apply the framework presented here to 
compute the baryon asymmetry in the two Higgs doublet model j^H] 
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